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ABSTRACT 

Ongoing star formation in the Milky Way requires continuous gaseous fuel from accretion. Previous 
work has suggested that the accretion of dwarf galaxies could provide the needed gas for this process. 
In this work we investigate whether dwarf galaxy accretion is consistent with the radial profile of star 
formation observed in the Milky Way, which is strongly concentrated toward the center of the gaseous 
disk of the Galaxy. Using realistic parameters for the Galactic potential, gaseous halo, Galactic disk, 
velocities of dwarf galaxies, and effects of drag on stripped gas in the halo, we model the distribution of 
radii at which dwarf galaxies accrete onto the Galactic disk. We also model the radial distribution of 
the accretion of gas that cools directly out of the halo by examining the results of recent simulations. 
We show that dwarf galaxies cannot reproduce the concentration of accretion toward the center of 
the Galaxy required by star formation. We also show that clouds that cool directly from the halo can 
reproduce this central concentration, and conclude that this difference is largely due to the discrepancy 
in absolute specific angular momentum between the two mechanisms. 

Subject headings: ISM: kinematics and dynamics, ISM: clouds, stars: formation, Galaxy: halo, evolu- 
tion 



1. INTRODUCTION 

Observations of young stars show that our Galaxy, like 
other spiral galaxies, is currently forming new stars out 
of the interstellar medium. It has been known for some 
time that this rate of star formation, > !Af Q /yr (Rana 



1991 hereafter Rana91), is so fast that it will exhaust all 
the gas in the star- forming area of the Milky Way fMW) , 



about 2x10 M(7) (Henderson et al 



(H 
!B 



1982 



Levine et al. 



1 200 6 1 hereafter LBHOb), in a small number of Gyrs. This 
timescale is much shorter than the age of the MW (> 10 
Gyrs), which leads us to conclude that we are either in a 
special moment in the evolution of the MW or that the 
MW is conti nuously gaining fuel fo r star formation. 

Work by |Chiappini et al.| (fl997| ) (hereafter CMG97), 
among others (see Larson 1972 for an early discussion), 
have bolstered the view that the MW is currently accret- 
ing new gas to form stars. These works showed that the 
distribution of metallicities of stars in the MW can be 
reproduced by invoking two epochs of gas infall: one in- 
tense, early epoch and one slower, ongoing epoch. These 
works also showed that star formation in the absence of 
the continuous accretion of fresh, relatively unenriched 
gas cannot generate the correct distribution of metals in 
the stars we observe: if all the gas had been available 
for star formation since z—l we encounter the "G-dwarf 
problem," wherein models predict too many long-lived 
stars with very low metallicities compared to observa- 
tions. 

In an extension of this type of s t ellar population 



metallicity modeling, Colavitti et al. (20081 (hereafter 



CMM08) showed that if gas accretion traces the accre- 
tion of smaller dark matter (DM) satellites (i.e., dwarf 
galaxies) in cosmologically simulated galaxies, the MW 
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distribution of metals in stars can also be reproduced. 
CMM08 also showed that this claim is only true in sim- 
ulated Galaxies with relatively quiescent recent merger 
histories, i.e., with all merger ratios below 1:20 in the 
last 8 Gyrs (since z = 1). These galaxies only accreted 
relatively small dwarfs during that time, consistent with 



the slower, ong oing accretion epoch from CMG97. Stew- 
lart et al. 1120071 have shown that about 40% of MW-sized 
halos in simulations of standard ACDM cosmology have 
merger histories that meet or exceed this recent merger 
quiescence criterion. This is to say that the constraints 
from CMM08 put the MW in the 40% most quiescent 
MW-sized haloes, and are therefore not necessarily un- 
reasonable. 

To do this modeling CMM08 had to assume that the 
accreted gas would be distributed within the Galaxy to 
match the observed radial distribution of star formation, 
as the DM-only cosmological simulations did not trace 
the evolution of the gas itself. By tracing the product s 
of star formation (Hn regions jG uestcn fc Mezgerf l 982 ) , 
Supernova remna nts ( |Leahy fc Wu||198H ), and pulsars 



(Lyne et al. 19851), it has been shown that star forma- 
tion in the MW is a strong function of Galactic radius, 
peaking at a Galactic radius R ~ 4 kpc, and tapering off 
to R ~ 15 kpc. The rate of this decrease varies depend- 
ing upon the star formation tracer employed, but it can 
be fit by an exponential function, 

T(R) = e^ R - R ®\ (1) 

where T is the star formation rate per unit area normal- 
ized to the solar circle (R = 8.5 kpc) and —0.35 kpc" 1 < 
< -0.18 kpc" 1 (see Figure [I}. 



It was suggested by |Maller~fc Bullock| ( |2004[ ) that a 
massive gaseous halo (MbaXo — 1O U -M ) around the MW 
could also be a source of fuel for ongoing accretion onto 
the Galactic disk. In t his theory, the hot (T ~ 10 6 K) 
halo unstably cools (see |Field||1965| |Murray fc Lin|1990 1 
to form overdense clouds that then fall onto the Galac- 
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Fig. 1. — The observed MW star formation rate per unit area 
normalized to the solar circle as measured by supernova remnants, 
Ho regions, and pulsars. Also shown are our fits to Equation [T] 
over the 4 kpc < R < 15 kpc range for each tracer. T his plot was 
adapted with permission from |Stahler fc Palla| [2005) ■ 



tic disk. Simulations have since confirmed that a mas- 



sive halo can generate such clou ds (e.g., Sommer-Larsen 
[2006] |Kaufmann et al.|2005| and |Kaufmann et al.|2 008a, 
hereafter KBMF08) and that this cooling halo material is 
consistent with the observed HVC s in the Galactic Halo, 
excepting the Magellanic Stream ( Peek et al.|2008 here- 
after PPS-L08). We therefore investigate these cooling 
halo clouds, alongside dwarf galaxies, as an alternative 
source of fuel for star formation in the MW. Indeed, it 
is not known which of these gas accretion mechanisms, 
dwarf galaxy stripping or the cooling of a massive halo, 
domi nates accretion onto spirals in the local universe 
(e.g., |Fraternali fc Binney||2008 1 . 

In this work we examine the accretion of gas from small 
dwarf galaxies and cooling halo clouds in terms of their 
consistency with the observed radial distribution of star 
formation in the MW. In particular, we make a two con- 
jectures about the radial movement of gas in the Galaxy 
(jp), which we use to compare the observed radial dis- 
tribution of star formation and the modeled radial dis- 
tribution of accretion. We construct a model of Galactic 
accretion, which is dependent upon a variety of physi- 
cal parameters of the Galaxy (O and the accreting gas 
(fQ. In fp] we discuss the relevant parameters of the 
Galax y, including the extent of gas in the Gal actic disk 
(j |3.1[ ), the gravitational potential of the MW (1 3.2 I and 
the structure of the Galactic gaseous halo ( |3 3[ ). In pi 
we discuss the kinematic and physical ch arac terist ics o f 
dwarf galaxies impinging on the MW (j ]4.1| and [4.2 I, 
meth ods by which gas can be stripped from these galax- 
ies ([4.3), and the characteristics of such gas once it is 
stripped (i 4.4 1. In fj5]we describe our method for model- 
ing gas accretion onto the Galactic disk. We discuss our 
findings in <|6]and conclude in f|7] 

2. RADIAL MOTION OF GAS WITHIN THE GALACTIC 

DISK 



As the MW presents a large area to accreting mate- 
rial (see [3.1 I, yet only forms stars towards the center, it 
is natural to ask whether gas accreted at one radius in 
the disk can propagate inward to form stars at another 
radius. Unfortunately, we are not provided with signifi- 
cant observational constraints. Given the radial surface 
density profile of gas in the MW (LBH06) , to fuel cen- 
tral star formation with gas from the Galaxy beyond the 
star-forming radius, an inflow velocity of order only 5 
km s -1 is needed. This fuelin g is too slow to de tect in 
observations of Hi in the MW. |Wong et al. (2004) inves- 
tigated inflow in other spiral galaxies, but were unable to 
put constraints on the inflow rate strong enough to rule 
out this relatively slow feeding. 

From a theoretical standpoint, to move material in- 
ward angular momentum must be redistributed. Indeed, 
it is energetically favorable for most of the disk gas to 
give its angular momentum to a small amount of mate- 
rial through some global mechanism. In this way most of 
the gas propagates inward while a small amount moves 
ever further from the center. The standard mechanism 
invoked for this angular momentum redistribution is spi 



ral density waves (e.g., Lynden-Bell & Kalnajs 1972 



hereafter L-BK72). To determine whether the torques 
generated by spiral density waves can generate signifi- 
cant inflow in the disk we define the inflow timescale, 
^inflow = L/t, where L is the angular momentum in 
the disk and r is the torque from spiral arms. We let 
L = TR 3 v c j > , and we use the d erivation of spiral torque 
from ( jBinney fc Tremaine 2008| ) for small pitch-angle spi- 
ral arms, 



ir 2 mRG(6^y 

P 



(2) 



where ST is the change in surface density in the spiral 
arms, k is the radial wavenumber (in rad -1 ) and m is the 
number of spiral arms. Inserting values for all of these 
parameters reasonable for the outer disk of the MW (see 
LBH06 for gaseous surface density maps) we find 



t= 1.2xl0 10 yr 



Rk 
2^ 



^ 4Af pc~ 2 



4M pc~ 2 J V ST 

approximately the age of the Universe. This is much 
longer than the timescale required for fueling in the disk, 
~ 10 9 years. 

Of course, we cannot conclusively claim that gas does 
not move radially from the outer disk to the areas of 
star formation. While torques from gaseous spiral arms 
do seem unlikely to be able to move gas on the relevant 
timescale, without empirical evidence we do not know 
whether some other mechanism might be at work. In- 
deed, other mechanisms have been proposed to induce 
significant radial inflow. We therefore consider two cri- 
teria, each of which could be true or false, regarding the 
movement of gas in the disk: the "weak no-flow crite- 
rion" and the "strong no-flow criterion." The weak no- 
flow criterion is that gas accreted beyond the maximum 
star- forming radius (R > 15 kpc) of the MW does not 
move into the star forming region. The strong no-flow 
criterion is that gas accreted at any radius in the MW 
outside 4 kpc forms stars only at the radius at which it 
lands. Since the MW bar, which could extend out to 4 
kpc, can be a very efficient mechanism for the redistri- 
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bution of angular momentum (e.g. Blitz fc Spergel|1991 
L-BK72), we make no assumptions about the region in 
terior to 4 kpc, and omit it from our analysis. 

We use these two criteria to test our models of Galactic 
accretion. If a model of MW accretion puts too much ma- 
terial in the outer disk, it will not meet the weak no-flow 
criterion. If a model of the MW accretion distributes gas 
along the disk with a surface density inconsistent with 
the observed profile of star formation, it will not meet the 
strong no-flow criterion. We note that the strong no-flow 
criterion is very stringent, and certainly not required by 
spiral torque theory. We include this as a criterion only 
to show whether models do or do not require movement 
within the star-forming disk. We describe the quantita- 
tive constraints on accretion that stem from these two 
criteria below. 

We know the amount of gas in the MW beyond the 
maximum star forming radius dominates the mass of 
stars in that region, M HItOUt ~ 1.7 x 10 9 M Q (LBH06), 
and that the current star-formation rate in the MW is 
> lM /yr. Since all accretion histories that are con- 
sistent with thin-disk stellar metallicities have a broadly 
declining ac cretion rate over the last 8 Gyrs (Chiappini 
et al. 19971, we can take a lower limit of the total star 
formation since z ~ 1 within the MW to be 8 x 10 9 Mq. 
A given ge neration of stars will return 48% of its gas 
to the ISM QRana||1991| ), which, after many generations, 
means that a unit of accretion will yield about 2 units of 
star formation. We therefore only require half the star 
formation rate to be represented by accretion onto the 
MW during this time, M acc = 4x 1O 9 M . If the amount 
of gas a model places beyond R = 15 kpc is too high 
compared with the amount it places within R = 15 kpc 
to match these observational constraints, the model is 
inconsistent with the weak no-flow criterion: gas must 
propagate inwards in the disk. We define this limiting 
fraction as a out = M HI , out /M acc = 0.44. 

Given the observations of star-formation tracers (Fig- 
ure [T| , we were able to parameterize the radial profile of 
star formation within the disk with Equation [l] For a 
model of accretion onto the MW to be consistent with 
the strong no- flow criterion, the radial profile of accre- 
tion in the model must fall within the range of these 
radial profiles of star formation. If at some place in the 
disk the accretion were to lag behind the star formation, 
the star formation would have ceased; if the star forma- 
tion were to lag behind the accretion we would see large 
gas reservoirs at that radius. Since neither of these is 
seen, any model which has a radial accretion profile with 
a (3 < —0.35 kpc -1 or (3 > —0.18 kpc -1 is inconsistent 
with the strong no-flow criterion. 

3. THE GALAXY 

3.1. The Scale of the Galactic Disk 

The extent of the Galaxy disk can be measured in 
terms of DM, stars, or gas. In this case we wish to de- 
termine the extent of the Galactic disk in terms of the 
ability to accrete (arrest) gas. Since inf ailing gas could 
pass by a disk of stars or of DM, but not of gas, it is 
the extent of the gaseous disk that is relevant to our disk 
size measurement. LBH06 show through observations of 
Hi that the MW disk has a surface density above 0.1 M 
pc -2 , or 3 xl0 19 cm -2 in hydrogen atom column density 



within Rd ~ 30 kpc. This is the relevant column density 
to arrest the infall of accreting material, as most accret- 
ing gas is expected to be near this column (see [4.4 1. If 
we choose instead a minimum surface density to arrest 
accreting gas of 10 20 cm -2 , we find a maximum radius of 
~ 25 kpc, which leaves the qualitative conclusions of this 
work unchanged. Therefore, in our modeling, we take the 
Galactic disk to be capable of accreting gaseous material 
out to R D = 30 kpc 

3.2. The Gravitational Potential of the Galaxy 

To understand the process of accretion in the MW one 
must first parameterize the dominant force acting on ac- 
creting material: gravity. The gravitational potential 
of the Galactic DM halo can be modeled by examining 
the movement of stars and gas in various positions in 
the MW and MW halo. We use the model from Wolfire 
al. (19951), consisting of a Miyamoto-Nagai disk po- 
tential, aspherical bulge potential and a logarithmic DM 
halo potential. These three potentials are added together 
and constitute the total Galactic potential we use in our 
models of accretion and the structure of the gaseous halo. 



3.3. The Cast 



Halo 



In this work, when we refer to the gaseous Galactic 
halo, we are referring to the gas that occupies the vol- 
ume surrounding the MW out to approximately the virial 
radius (R V i r ~ 250 kpc), where it is assumed to merge 
with the inter-galactic medium. It is this halo gas that 
will strip gas out of dwarf galaxies and generate drag on 
these stripped clouds. It is therefore crucial to determine 
the extent and character of the gaseous halo. 



Observations th e soft X-ray background (e.g., Kuntz 
& Snowden 20001, the pres sure confinement of high- 
velocity clouds (HVCs; e.g., Wolfire et al. |1995|, and 
absor ption of high ions in X-ray spectra (e.g., |Fang et al. 
2006) suggest the the existence of such a halo. Gon- 
stramts on the X-ray emission measure from the halo 
(4.7 xlO -3 cnr 6 pc) and temperature constraints from 
both the soft X-rays and absorption by high ions (10 6 K 
< Thaio < 3 x 10 6 ) allow us to construct a simple hy- 
drostatic, isothermal halo within the Galactic potential 
that depends upon only the temperature of the halo. Un- 
fortunately, this relatively narrow range in temperature 
can significantly affect the mass of the halo and density 
of the halo beyond ~ 30 kpc. We therefore use one of two 
distinct models of the halo in each of our runs: a "low- 
mass" halo, with T = 10 6 K and M haio ~ 1O 9 M , and a 
"massive" halo with T=3 x 10 6 K and M halo ~ 1O U M . 
We note that while this is certainly a simplistic formu- 
lation of the halo, the massive halo is consistent with 
the density structure of the simulated halo described in 
Kaufmann et al. (2008b I to within a factor of 2. 



4. DWARF GALAXIES 

4.1. Kinematic Characteristics of Inf ailing Dwarf 
Galaxies 

To determine the distribution of gas deposition on the 
Galactic disk by dwarf galaxies we must first determine 
their initial conditions as they fall into the Galactic po- 
tential. We would prefer to determine this by observa- 
tional evidence, but the observed dwarfs around the MW 
may not be a good representation of typical dwarfs enter- 
ing the influence the MW for the first time. Dwarfs on a 
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"collision course" with the core of the MW disk may only 
be observed once near the virial radius, whereas dwarfs 
that avoid such collisions may be observed on many sub- 
sequent orbits. More generally, the ellipticity and veloc- 
ity of dwarf galaxy orbits will determine the amount of 
time a given galaxy spends near the virial radius. Addi- 
tionally, we do not have a f ull census of dwarf g alaxies 
at the MW virial radius (see Koposov et al.|20"07 l. 

Instead of using observational data, we use radial and 
azimuthal v elocity distribu tions determined from DM 
simulations. Benson (2005) compiled data from a num- 
ber of different DM cosmological simulations to deter- 
mine the distribution of radial and tangential veloci- 
ties of satellites joining larger halos at the virial radius, 
/ (v r , v^). He found that / (v r , v^) is relatively constant 
over a wide range of primary halo masses and cosmic 
time. We draw our sample of dwarf galaxy velocities 
from this probability distribution function to test our 
accretion hypotheses. We note that this distribution of 
subhalo velocities applies both to the dwarf galaxies we 
wish to examine and and DM haloes that may not con- 
tain galaxies. 

4.2. Mass of In] 'ailing Dwarf Galaxies 

Dwarf galaxies around the MW and the local group 
have a large range of sizes, starting with the Large Mag- 
ellanic Cloud (LMC) with M ~ 2 x 1O 1O M and going 
beyond our observational limits a t near 10 6 M Q (Koposov 
et al.|2007[ |S mon fc Geha|2007| . Since our observations 
are rather incomplete, we cannot determine the relative 
importance of large and small satellites in fueling accre- 
tion onto the MW. While both extended Press-Schecter 
theory and simulations show that, per decade, massive 
satellites contribute most to the accretion rate of MW- 



sized Galaxies (e.g., Fakhouri & Ma 2008 Stewart et al. 
2007 and references therein), they also show that much of 
the accretion is generated by the ma ny decades of smaller 
mass satellites. Madau et al. ( 2008 1 show that sub-halos 
in very high resolution DM simulations have equal DM 
mass per decade, but they do not address the question 
of halo accretion directly. Given this uncertain situa- 
tion we cannot limit our investigation to only massive 
or low-mass satellites, as either may be the dominant 
source of Galactic fuel. To test this broad range of dwarf 
galaxy masses, we run each model with a single dwarf 
galaxy mass parameter, which we allow to vary over a 
large range from model to model. As we will see in the 
next section, dwarf galaxy mass determines when the gas 
becomes stripped from the stars. 

4.3. Gas Stripping and Drag in Dwarf Galaxies 

Until gas is stripped out of dwarf galaxies it has little 
chance of entering the disk of the Galaxy, as dwarf galax- 
ies will retain the majority of their angular momentum, 
and typically have perigalactions > 30 kpc. Gas can 
be removed from galaxies (without being destroyed) by 
either ram-pressure stripping or tidal stripping. In ram- 
pressure stripping, the pressure exerted by the halo gas 
"head-wind" on the dwarf galaxy gas exceeds the grav- 
itation force per unit area on the gas from the dwarf 
galaxy itself. In tidal stripping, the differential force on 
two ends of the dwarf galaxy from the MW gravity tears 
the galaxy apart. This is facilitated by generally larger 
size of the gaseous component of gas-rich dwarf galaxies 



as compared to the stellar component. In the case of the 
mergers of interest, we will neglect tidal stripping as a 
method for removing gas from dwarf galaxies and focus 
only on the effects of ram-pressure stripping. We will 
verify the inconscquentiality of this assumption on our 
results in Sj6] 

It is important to realize that some or all of this 
stripped gas could eventually shock, become hot and 
merge with the gaseous halo. We cannot distinguish in 
our modeling between (a) gas that was originally part 
of a dwarf galaxy, merged with the halo, and later 're- 
cooled' via the process described in <|l] and (b) gas that 
cools from the halo having never been a part of a dwarf 
galaxy. Given this, we allow the dwarf galaxy accretion 
scenario the 'benefit of the doubt' and assume that ma- 
terial does not dissipate into the Galactic halo. This 
allows us to ask whether dwarf galaxies can fuel Galactic 
accretion in the absence of a massive, cooling halo. 

To determine when a dwarf gal axy becomes s t ripped 
we use the criterion established in Gunn & Gott (19721. 
Their stripping criterion can be formulated as 



c 



1 



"dwarf C 2 < "-halo" 2 , 



(4) 



where "halo is the number density of the halo at the lo- 
cation of the dwarf, v is the velocity of the dwarf with 
respect to the halo gas, a is the velocity dispersion of the 
stars in the dwarf galaxy, and "dwarf is t he typical num 
ber density of gas in the dwar f gal axy (see Grevich & Put 
manj2008 ) . As we stated in {A.2 there is a large range of 



possibly relevant dwarf galaxy mass, so we choose three 
scenarios to test: one in which most satellites are so small 
that they immediately exceed the criterion in Equation 
[4] upon entering the halo (( = (km s -1 ) 2 cm -3 ), one 
m which most dwarf galaxies have a — 10 km s -1 and 
"dwarf = 1 cm -3 (£ = 33 (km s -1 ) 2 cm™ 3 ), and one in 
which the dwarf galaxies never become stripped (£ = oo 
(km s -1 ) 2 cm -3 ). These extremal cases are not intended 
to represent reality, but rather act as tests of the sensi- 
tivity of the model results to the nature of gas stripping. 

Once gas is stripped from the dwarf, ram pressure drag 
from the gaseous halo of the MW can become a signif- 
icant force on the resulting cloud. To model the effect 
of this drag on the cloud we use the simple ram-pressure 
drag expression, 



1 



rPhalo-u C d A, 



(5) 



where A is the ar ea of the cloud and Cd is th e drag coef- 
ficient. Following Benjamin & Danly (1997), we assume 
Cd = 1. This equation is equivalent to 



ad 



(6) 



2N H ' 

where a d is the acceleration of the cloud and Njj is the 
hydrogen column number density of the cloud. 

4.4. The Structure of Stripped Gas 

In the Galactic halo we have evidence of stripped gas 
from dwarf galaxies in observations of the Magellanic 
Stream. The Magellanic Stream is a long chain of gas, 
detected primarily in Hi as HVCs. HVCs are seen all 
over the sky, but the Magellanic Stream is the only 
undoubtable case of dwarf galaxy stripping around our 
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TABLE 1 
Parameter list 



19 20 21 

logio ( N Hi) [logio (cm" 2 )] 

Fig. 2. — The distribution of mass of the Magellanic Stream in 
Hi by logarithmic column density, assumi ng that the Stream lie s 
70 kpc from the observer. Data are from |Putman et al.| ( |2003b| l . 
The mass-weighted column density of the Stream peaks nearlfj 30 



Galaxy. Therefore we we use the Magellanic Stream to 
determine the typical column density of material that 
has been stripped from dwarf galaxies. 

Unfortunately, there seems to be no specific preferred 
column density in surveys of the Magellanic Stream or 
HVCs in general (Wakker 2004). HVC columns begin 



above 10 20 cm" 2 , but have a lower bound that is de- 
termined by the sensitivity of the observations. For the 
question of accretion, though, we are not concerned with 
the typical column density, but rather the typical mass- 
weighted column density: 



EZi N HA N Hti ilD 2 fj,m p 



(Nh) 



(7) 



where i ranges over resolution elements, M is the num- 
ber of such elements, Nhs is the column density of each 
element, m, is the mass of the cloud at each element, f2 
is the solid angle subtended by such an element, and D 
is the distance to the cloud. This mass- weighting effec- 
tively diminishes the significance of the long tail towards 
low column density, which contributes little to the overall 
accretion rate. For the observations o f the entire Magel- 
lanic Stream by Putman et al. (2003b) at 15'. 5 resolution, 
(Nh) = 10 19 6 cm _i , with most of the mass of the Stream 
residing between 10 194 cm -2 and 10 20 cm~ 2 (see Figure 

We note that due to resolution effects the measured 
(Nh) is a strict lower limit. This is to say that it is 
possible that the Magellanic Stream is made up of many, 
higher column density clouds unresolved by observations, 
which would yield a higher true (Nh), although HVC 
absorption- line observations of distant quasars show that 
the column d ensities are relative ly converged at a reso- 
lution of 10' ( fWakker et aT1|200T] ). 

We may also want to consider gas accretion from DM 
subhaloes that have not formed stars (dark galaxies). 
Modern simulations predict that ^70% of the DM ac- 
creted in subhaloes will have masses above 10 6 M M and 



parameter 


value / reference 


Rd 


30 kpc 


Rvir 


250 kpc 




IWolfire et al.| (11995 ) 


/ (v r ,V^) 


Benson[^ote|| 


N H 


> TO 19 '* cm - * 


M halo 


10 9 iW Q , 10 n M o 


^halo 


10 6 K, 3 x 10 6 K 


c d 


1 



Note. — A summarized list of parameters used in our models 
or the reference to the functions used. 

that there are of order 15000 such h aloes surrounding the 



MW , far more than are o bserved (Springel et al. 2008 



e. g. Simon fc GehapOOT I . Were these haloes to contain 



condensed gas, but not stars, they could be a dominant 
accretion mechanism of the Galaxy, and, due to their 
size and starlessness, could have significantly different 
gas distribu t ions t han are observed in the MS. [West- ' 



meier et al. (20051 conducted observations of outskirts 
of Andromeda and ruled out gas in such dark galaxies 
above 6 x 1O 4 M in a region that should contain dozens 
of galaxies with Mdm > 1O 6 M and many galaxies with 
Mdm > 1O 7 M . This puts the gas fractions of possible 
dark galaxies below 10~ 3 , f ar below the typical ~ 10 _1 
of standard dwarf galaxies (Mateo 1998), implying that 
if such subhaloes do exist they do not contribute signif- 
icantly to accretion of condensed gas onto the MW. By 
contrast, the lowest mass known gas-bearing galaxy, Leo 
T, contains significant dense Hi, with most mass resid- 
ing at column densities above 10 20 cm~ 2 . We therefore 
do not further consider the effects of the accretion of gas 
via dark galaxies. 

This discussion has only addressed the question of 
stripped neutral gas, mostly because observations of dif- 
fuse Hi are much easier to conduct that of diffuse Hn. 
Observations in H a of the Magellanic Stream by |Put-| 



man et al. (2003a 



have shown a significant amount of 
warm Hn gas exists in this system, although estimating 
the mass or column density of warm Hn from these mea- 
surements is diffi c ult. S imulations conducted by [Bland- 



Hawthorn et 
the Stream 



al. 



(2007) reproduce the observed Ha in 
Tneir results show that the typical mass- 
weighted column density of warm Hn is above 10 20 cm" 2 
and exceeds that of the Hi. We therefore tentatively con- 
clude that in the case of gas stripped from dwarf galax- 
ies, the contribution from warm Hn can only increase the 
typical column density from the value measured from Hi, 
and we therefore consider our Nh estimate from the Hi in 
the Magellanic Stream to be a lower limit. 

5. METHODS 

Given these parameters of the disk, halo, and satellites 
of the MW (summarized in Table [I]), we construct a first- 
order finite-difference code to model the orbits of dwarf 
galaxies and clouds within the halo. The code evolves the 
position and velocity of a large number of objects, each 
representing a cloud or a dwarf galaxy around the MW. 
Galaxy objects are subject to only the force of gravity, 
V4>, from the MW. Cloud objects are subject to both 
gravity and ram pressure drag (as in Equation [6]) . Note 
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that this code does not capture any of the more com- 
plex physical processes that are modeled in modern hyr- 
drodynamical simulations. We believe that our simpli- 
fied treatment is an appropriate model for this unstudied 
question. 

In our model runs of dwarf galaxy accretion, 10 4 galaxy 
objects are started at random positions at the virial ra- 
dius wi t h velo cities drawn from the distribution from 
Benson| ([2005). These galaxy objects orbit until the 
stripping criterion is met (Equation W| , at which point 
we convert the galaxy object into a cloud object. We 
follow the orbits of these cloud and galaxy objects until 
they cross the disk of the galaxy within 30 kpc, or until 
8 Gyrs pass in the the model. 

In our model runs of condensing cloud accretion we use 
the velocities and positions from clouds found in simu- 
lations from PPS-L08 and KBMF08 as our initial condi- 
tions, and evolve them as cloud objects as we would in 
the dwarf galaxy accretion model. Note that we do not 
use any information from these simulations other than 
the positions and velocities of the clouds selected at an 
arbitrary time-slice (0 < z < 1) to give a comparable 
treatment to the dwarf galaxy model. 

The results of the model runs are radii at which each 
object lands in the the Galactic disk. Since the to- 
tal number objects and masses of the objects are ir- 
relevant to the physics of the model, we do not gen- 
erate any constraints related to the overall accretion 
rate of the disk. The overall accretion rate is a com- 
pletely free parameter in our model; readers interested 
in the consistency between the MW star-formation rate 
and the overall accretion rate of cooling halo clouds and 
dwarf galaxies should refer to PPS-L08 and GP09, re- 
spectively. Instead, the primary results of these model 
runs are the relative accretion rates as a function of ra- 
dius. We use this accretion as a function of radius to 
test whether each model meets the weak no-flow crite- 
rion (a ou t < 0.44) and whether it meets the strong no- 
flow criterion (-0.35 kpc" 1 < (3 < -0.18 kpc" 1 ). 

Given the uncertainties in the hal o pr ofile ( |3.3| , the 
relevant cloud co lumn density Nh ( §4.4fl , and the strip- 
ping criterion £ ({ 4.3 ), we have a variety of runs for each 
combination of these criteria. To simplify the interpre- 
tation of the results, for each run we choose single values 
for N H (10 19 - 6 or 10 20 6 cm" 2 ) and C (0, 33, or oo (km 
s -1 ) 2 cm -3 ), which we apply to all clouds and galaxies in 
a given run. Results are shown in Tableland discussed 
below. 

6. RESULTS 

In Table [2] we show the results from the dwarf galaxy 
and cooling cloud accretion models. For each model 
we show two parameters. The first parameter, a out , is 
the ratio of material accreted beyond 15 kpc to that of 
the material accreted within 15 kpc, which determines 
whether a model meets the weak no-flow criterion. The 
second parameter is the coefficient of the best exponen- 
tial fit to the accretion rate radial profile between 4 and 
12 kpc ((3, as in Equation [lj , which determines whether 
a model meets the strong no-flow criterion. All 12 dwarf 
galaxy accretion models fail to meet either the weak or the 
strong no-flow criteria. It is also interesting to note that 
only in the simulations with £ = (km s" 1 ) 2 cm~ 3 and 
Nh = 10 19 ' 6 cm -2 did all of the gas accrete onto the disk 



within 8 Gyrs. This result is consistent with the supposi- 
tion that stripping and drag are crucial to gas accretion 
from satellites. We also note here that adding further 
stripping through tidal effects cannot change this result, 
as even using immediate stripping the modeled accretion 
does not meet the criteria. 

By stark contrast, all of the models generated from 
simulated data of cooling haloes meet both the weak 
and the strong no-flow criteria. In a cooling-halo sce- 
nario, therefore, no radial gas movement within the disk 
is needed to reproduce the observed star formation ra- 
dial profile. Given the range in the observational data, 
we cannot choose one of the models as a better fit, and we 
note the surprising similarities between the results from 
the PPS-L08 data set and KBMF08 data set. These 
simulations were run entirely independently by different 
investigators with unrelated initial conditions and differ- 
ent resolutions, and yet they produce very similar results 
in terms of their accretion profile in our models. 

In Figures [3] and [4] we show the results of a larger num- 
ber of dwarf galaxy model runs. Figure [3] shows whether 
models meet the weak no-flow criterion (a ou t < 0.44), 
Figure [4] shows whether models meet the strong no-flow 
criterion (-0.35 kpc" 1 < f3 < -0.18 kpc" 1 ). The plots 
show quite clearly that, for both haloes, only runs that 
have unreasonably small values of Nh and ( < 20 (km 
s" 1 ) 2 cm" 3 can meet either of the no-flow criteria. 

In Figure [5] we show the equivalent results for the cool- 
ing halo models. Here, the stripping criterion, £, is ir- 
relevant (cooling clouds are not associated with stars or 
DM haloes), so we plot both (3 (left) and a out (right) for 
each model, using both the PPS-L08 and KBFM08 sim- 
ulations. We find that both the PPS-L08 and KBFM08 
simulations give consistent results in a ou t for all values 
of Nh and reasonable results for j3 for many values of 
Nh- For unreasonably low values of Njj, the cooling 
halo clouds tend to have too steep an accretion profile 
{13 < -0.35). 

The sharply differing results between dwarfs and 
cooled halo clouds are not surprising when considering 
the absolute specific angular momentum, |Z|, of the ini- 
tial objects. In both the PPS-L08 and KBFM08 simula- 
tions the (\l\) of the clouds is ~ 2.5 x 10 3 km s" 1 kpc, 



whereas the dwarf galaxies from Benson ( 2005 1 had 
(|/|) ~ 2.5 x 10 4 km s _1 kpc. The gas in dwarf galaxies 
had vastly more |/| to get rid of before reaching the disk, 
much increasing the chances of hitting the outside of the 
disk before reaching the center. 

7. CONCLUSIONS 

We conclude from this work that the accretion of gas 
from dwarf galaxies cannot be the primary source of fuel 
for star formation in the MW without significant inward 
movement of gas throughout the Galactic plane. We also 
conclude that simulations of material cooling out of the 
Galactic halo naturally generate a population of clouds 
that is consistent with the radial profile of star formation 
from R = 4 to R = 15 kpc. We find that a natural 
distinction between these two accretion mechanisms is 
their typical specific angular momentum. 

We hope this work will spur investigation into the fun- 
damental parameters of ongoing Galactic accretion. One 
important question this work raises is whether halo drag 
on accreting cold gas can be modeled with the simple 
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18.0 18.5 19.0 19.5 20.0 18.0 18.5 19.0 19.5 20.0 

kgio ( N m ) [iogio (cm" 2 )] log 10 (N m ) [log 10 (cm" 2 )] 

Fig. 3. — Two plots of a large number of dwarf galaxy accretion model runs, the left plot representing runs with a low-mass halo, the 
right plot representing runs with a high-mass halo. Each pixel represents a run with 10 4 dwarf galaxy objects; we vary the the column 
density of accreting gas (Nh) on the x-axis and the stripping criterion of the gas on the y-axis. Dark pixels represent model runs that 
meet the weak no-flow criterion (a ou t < 0.44). We examine a much larger range in Nh than we consider observably reasonable in or der to 
determine our sensitivity to this assumption. We also show our lower limit on typical Nh with the vertical line in each plot (see | ]4.4| l. This 
plot demonstrates that for dwarf galaxy accretion to meet the weak no-flow criterion, typical gas column density needs to be unreasonably 
low. 
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Fig. 4. — Two plots of a large number of dwarf galaxy accretion model runs, the left plot representing runs with a low-mass halo, the 
right plot representing runs with a high-mass halo. Each pixel represents a model run with 10 4 dwarf galaxy objects; we vary the the 
column density of accreting gas (Nh) on the x-axis and the stripping criterion of the gas on the y-axis. Dark pixels represent models that 
meet the strong no-flow criterion (—0.35 kpc -1 < f3 < —0.18 kpc — 1 ). We examine a much larger range in Njj than we consider observably 
reasonable in order to det ermine our sensitivity to this assumptions. We also show our expected lower limit on typical Njj with the vertical 
line in each plot (see j ]4.4| . This plot demonstrates that for dwarf galaxy accretion to meet the strong no-flow criterion, typical gas column 
density needs to be unreasonably low. 
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TABLE 2 

Model run results 



Run# 


method 


halo mass 


c 

(km s 1 ) cm 3 


log iVff, c loud 

log cm - 2 


Clout 


& 

1/ kpc 


/accr 


I 


Hwa VTCi 


1UW IIIOjOD 




19.6 


2.9 


_n m 4- n 02 


0.09 


2 


dwarfs 


low-mass 


oo 


20.6 


3.1 


0.02 ±0.02 


0.09 


3 


dwarfs 


low-mass 


33 


19.6 


2.7 


-0.01 ± 0.03 


0.09 


4 


dwarfs 


low-mass 


33 


20.6 


2.9 


0.05 ±0.03 


0.08 


5 


dwarfs 


low-mass 





19.6 


3.9 


0.10 ±0.01 


1.00 


6 


dwarfs 


low-mass 





20.6 


4.0 


-0.00 ±0.02 


0.19 


7 


dwarfs 


massive 


oo 


19.6 


2.5 


0.02 ±0.02 


0.09 


8 


dwarfs 


massive 


oo 


20.6 


2.8 


-0.02 ± 0.02 


0.09 


9 


dwarfs 


massive 


33 


19.6 


3.0 


-0.07 ±0.02 


0.09 


10 


dwarfs 


massive 


33 


20.6 


2.8 


-0.01 ± 0.02 


0.09 


11 


dwarfs 


massive 





19.6 


4.5 


0.10 ±0.01 


1.00 


12 


dwarfs 


massive 





20.6 


4.5 


-0.02 ±0.02 


0.19 


13 


PPS-L08 


massive 


N/A 


19.6 


0.07 


-0.20 ±0.02 


1.00 


14 


PPS-L08 


massive 


N/A 


20.6 


0.19 


-0.13 ±0.03 


1.00 


15 


KBMF08 


massive 


N/A 


19.6 


0.13 


-0.27 ± 0.02 


1.00 


16 


KBMF08 


massive 


N/A 


20.6 


0.26 


-0.27 ±0.02 


1.00 



Note. — The input parameters for each of the models we run to the left of the divider, with the results at the right of the divider. 
The large number of models is require d to vary each model parameter independently. The method column denotes whether we generate 
the clouds from dwarf galaxies (as in j ]4.1[ l or from one of the two cooling halo simulations. The £ column denotes at what point gas is 
stripped from dwarf galaxies, as described in j ]4.3| and Equation |4] The log Nu jC i ou d column denotes the assumed column density of the 
clouds. We measure the steepness of the normalized accretion rate per unit area with radius (/3, see Equation!]! and the ratio of material 
accreted beyond 15 kpc to that accreted within, a out . We also measure the fraction of clouds or galaxies accreted onto the disk within 8 

Gyrs, /accr- 




18.0 18.5 19.0 19.5 20.0 18.0 18.5 19.0 19.5 20.0 

logio (Nm ) [log 10 (cm" 2 )] log 10 (N m ) [bg 10 (cm" 2 )] 



Fig. 5. — The two plots are of /3 at left and a ou t at right as a function of gas column density (Nh) for the data taken from the cooling 
halo simulations (PPS-L08 in gray, KBMF08 in black). The gray box in the left plot represents the range in j3 found in star-formation 
tracers (the strong stay put criterion), the line at the top of the right plot shows the maximum allowable value of a out consistent with the 
weak no-flowno-flow criterion. The lower limit for Nh determined in §4.4| is also shown as a line at the right of each plot. 
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ram-pressure drag formulation expressed in Equation [6] 
The complex interface between l ow- and high-densit y gas 
has bedeviled simulations (e.g., Agertz et al.||2006 l, and 
may produce different drag effects than we have modeled. 
If the drag were significantly stronger our first conclusion 
could be incorrect, although modern grid-based simula- 
tions seem to find d < 1 (Heitsch et al, in prep). A 
concerted effort to determine whether ram-pressure drag 
is a relatively good assumption for neutral clouds in the 
halo and what an appropriate choice for Cd is in this 
context could al so be of use to studies of HVCs (e.g., 
Peek et aT||2007 1 . Another question raised by this work 
is whether all cooling halo models generate results con- 
sistent with those found from PPS-L08 and KBMF08. 
While no parameters were tuned in either of these sim- 
ulations to match the star formation data, we cannot 
claim conclusively that all cooling halo models will re- 



produce these results, and further investigation by other 
groups is needed. Along these same lines, the overall 
accretion rate from these simulations must be examined 
in more detail. The accretion rate in PPS-L08 is only 
~ 0.2Af Q /yr, significantly below that of the observed 
MW rate. It is currently unknown whether that result 
is a robust prediction of all cooling-halo models, but if 
it were confirmed it would cast doubt on the ability of a 
cooling halo scenario to fuel MW star formation. 

The author would like to thank Tobias Kaufmann and 
Jesper Sommer-Larsen for supplying simulation data for 
the cooling halo models, Eugene Chiang for the calcula- 
tion in ^2] and Mary Putman, Karin Sandstrom, Steve 
Stahlcr, Louis-Benoit Desroches and Kathryn Peek for 
helpful comments. The author was supported in part by 
The National Science Foundation grants AST-0709347 
and AST-0406987. 
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